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Abstract 

We study a large sample of RXTE PCA/HEXTE observations of 
Cyg X-l. We characterize the spectra by soft and hard X-ray colours 
(which define the spectral states), and fit them with a physical model of 
hybrid, thermal/non-thermal Comptonization. We then fit the power 
spectra by a sum of Lorentzians. We show the resulting correlations 
between the spectral colour (defining the hard, intermediate and soft 
spectral states), the Comptonization amplification factor, the frac- 
tional Compton reflection strength, and the peak frequency of the 
Lorentzians. We also calculate the fractional variability (rms) as a 
function of photon energy. We fit the obtained rms dependencies by 
physical models of varying accretion rate, a varying soft-photon input 
from an outer disc to a hot inner plasma, and a varying dissipation 
rate in a hot corona above a disc. 



1 Introduction 

Cyg X-l is an archetypical high- mass X-ray binary containing a persistent 
black-hole source accreting from an OB supergiant. See, e.g., Ziolkowski 
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Figure 1: (a) The colour-colour diagram for our sample (black points; see Section 
[1] for colour definition) . The blue line specifies our definition of the spectral state 
parameter, changing from to 1. (b) The histogram of the observation number 
(normalized to unity) as a function of the state parameter. 

(2005) for the properties of the binary, and, e.g., Zdziarski & Gierlinski 

(2004) for a review of its spectral and timing properties. 

This work is, in some part, an application to Cyg X-l of the timing 
and spectral analysis methods developed for the black-hole low-mass X- 
ray binaries XTE J1550-564 and XTE J1650-500 by Gierlinski & Zdziarski 

(2005) . An independent analysis of the spectral correlations only in Cyg 
X-l is that of Ibragimov et al. (2005). In addition, we fit the power spectra 
of Cyg X-l by a sum of Lorentzian components and look for correlations of 
the frequencies with other parameters. See Axelsson, Borgonovo &: Larsson 
(2005) for another study of that type. 

We can distinguish two approaches to studying X-ray data. In the first 
one, one may study a small selected sample in great detail. The advantage 
of this approach is the possibility of studying physics of the sources in detail, 
and to determine their emission mechanisms and geometry. A disadvantage 
is that it gives little information regarding the source variability, or dynam- 
ics. In the second approach, one may study a large sample of observations 
using approximate methods, e.g., use a colour-colour diagram. Although 
this is rather phenomenological, it allows us to study the source behaviour 
over all its states and it gives some indications about the source dynamics. 

In the present analysis, we combine these two approaches. We define here 
a representative sub-sample chosen from the large number > 700 pointed 
observations of Cyg X-l by RXTE performed up to 2004. The chosen sub- 
sample covers all spectral colours and all types of variability. Then, we fit 
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Figure 2: (a) The colour-colour diagram of all observations (red points) showing 
the subsample studied in detail (33 black circles), (b) The fractional variability as 
a function of the state parameter. 



the selected data in detail by a physical model. We show the data form a 
well-defined sequence as a function of the spectral state. This allows us to 
show a number of clear correlations between various fitted parameters. 

In Fig. Ufa), we show the positions of all of the considered observations 
in a hard-soft colour diagram, based on our fits to the RXTE PCA data. 
The colours are defined as in Done & Gierlihski (2003). The soft and hard 
colours are defined as the ratios of the intrinsic, absorption-corrected, energy 
fluxes in the energy bands 4-6.4 keV to 3-4 keV, and 9.7-16 keV to 6.4-9.7 
keV, respectively. We see that the observations form a straight line on this 
plot, with relatively little scatter. This allows us to define our spectral-state 
parameter as the position on this diagram, changing from at the hardest 
spectra to 1 at the softest spectra. Fig. QJb) shows the histogram of the ob- 
servation number vs. the state parameter. We see two peaks corresponding 
to the hard and soft states. 

Then we define our smaller representative sub-sample, containing 33 
data sets studied in detail. They are marked on the colour-colour diagram 
by black circles in Fig. 12(a), and also in Fig. [2(b), which shows the fractional 
rms variability of the PCA light curves as a function of the state parameter. 
We see the rms first declines from the hard state toward softer ones, but it 
increases again in the softest states. See, e.g., Zdziarski & Gierlihski (2004) 
and Done, Gierlihski & Kubota (2007) for discussions of properties of the 
spectral states. 

Our spectral model is that of hybrid Comptonization, eqpair (Coppi 
1999; Gierlihski et al. 1999) in xspec (Dorman & Arnaud 2001). The 
hybrid electron distribution in this model consists of a Maxwellian and 




Figure 3: The six characteristic spectra of Cyg X-l (in the dF/d In E = EFe repre- 
sentation), from (a) very hard to (f) very soft, fitted by the hybrid Comptonization 
model. See Section [2] for details. 



a non-thermal tail. The model is characterized by the two main param- 
eters, the soft (disc) luminosity (or compactness), £ s , and the hard (hot 
plasma) luminosity/compactness, £^. The compactness is defined as usual, 
I ~ Lai/m e (? : where L is the luminosity, <tt is the Thomson cross sec- 
tion and m e is the electron mass. Hereafter, as our main parameter we use 
the Comptonization amplification factor, i.e., the hard-to-soft ratio, £h/£s, 
which is closely related to the X-ray (photon) spectral index, T. Then, the 
hard and intermediate states correspond to £^/£ s <^ 1, and the soft state 
corresponds to £h/£ s 1. 



2 Results 

Figs. Ela-f) present six selected representative energy spectra, EFe, in a 
sequence from the hard to the soft ones. The numbers in each panel give 
the RXTE observation ID. The red curves show the un-scattered blackbody 
component. Scattered blackbody photons provide then seeds for Comp- 
tonization by the hybrid plasma. The Comptonized component is shown by 
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Figure 4: (a) The Comptonization amplification factor, £h/£s, vs. the spectral 
state. The black filled squares and magenta empty squares correspond to the hard 
and intermediate state, respectively. The blue filled triangles and the red empty 
squares correspond to the soft state with a low, < 0.22, and high, > 0.22 values of 
the total rms. (b) The frequencies (in Hz) of the Lorentzian components, fitted 
to the power spectra vs. the spectral state. 



the blue curves. Then, that emission is Compton reflected from an optically- 
thick accretion disc (Magdziarz & Zdziarski 1995), which also gives rise to 
the Fe K line emission, the sums of both are shown by the green curves. 

Fig. Ufa) shows the dependence of the fitted Comptonization amplifica- 
tion factor, £h/£ s , on the state parameter. We see a very clear anticorre- 
lation, with our hardest spectrum having £h/£ s — 10 and the softest one, 
£h/£s — 0.25. Here, we define the hard state by £h/£s > 3.5, the intermediate 
state by 1 < £h/£s < 3.5, and the soft state by £h/£s < 1- In addition, we 
distinguish the soft states with low and high variability, see Fig. [10] below. 

Fig. [5ja) shows the fractional strength of the Compton reflection (which 
can be approximately identified with the solid angle subtended by the re- 
flection normalized to 2tt, R = Q/2ir) as a function of the Comptonization 
amplification factor. We see that the reflection strength decreases with the 
increasing spectral hardness, which is in agreement with previous results 
(Ueda, Ebisawa & Done 1994; Zdziarski, Lubinski & Smith 1999; Gilfanov, 
Churazov & Revnivtsev 1999, 2000; Revnivtsev, Gilfanov & Churazov 2001; 
Zdziarski et al. 2003). Fig. E^b) shows that the ionization parameter of the 
reflector decreases with the increasing spectral hardness, which is consistent 
with both the luminosity and the number of soft photons decreasing with 
the hardness. 

Figs. Ela-f) show the power spectra of the same selected 6 observations 
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Figure 5: (a) The relative strength of Compton reflection and the ionization param- 
eter of the reflector vs. the amplification factor, (b) The ionization parameter 
of the reflecting medium, £ = A~KF- lcm /n, where F- lon is the irradiating ionizing flux 
and n is the density of the reflector. 




Figure 6: The power spectra (in the dP/dln/ = fPf representation) for the six 
selected observations for which the energy spectra are shown in Fig. [3l The power 
spectra are fitted by up to four Lorentzians and an e-folded J -1 component. 
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Figure 7: (a) The frequency (in Hz) of the first Lorentzian, /i, vs. the Comp- 
tonization amplification factor, (b) The reflection scaling factor vs. f\. 



as those shown in Fig. [3l Each spectrum is fitted by a number of Lorentzian 
components and, in addition, an e-folded / _1 component, accounting for the 
power at low frequencies. The relative importance of the e-folded component 
increases with the increasing spectral softness. In Fig. |4^b), we see that all 
the Lorentzian components have the peak frequencies increasing with our 
spectral-state parameter, i.e., they are the highest for the softest spectra. 

We now show in Fig. EKa) the frequency of the first Lorentzian (/i, red 
in Fig. and in Fig. [6]) as a function of the Comptonization amplification 
factor, i\y/i s . We find a very clear correlation of /i oc (4/4)~ 2/3 . Fig. E^b) 
shows the relationship between fi and the strength of reflection, R. We see 
R increasing with fx, which is expected, given that R increases with the 
spectral softness (e.g., Zdziarski et al. 1999). 

Fig. Eta) shows the total rms (obtained by integration over the power 
spectra) as a function of the amplification parameter (or, equivalently, the 
spectral hardness). We see the rms is first decreasing as going from the 
hard (on the right) to the intermediate state, but then it increases again in 
the softest states. Here, we define the high-rms (> 0.22) soft state and the 
low-rms (< 0.22) one. 

Fig. 0(b) shows the total rms vs. the reflection strength. We see the 
reflection strength is increasing with decreasing rms, as the system moves 
from the hard to the soft state. However, the high-rms soft state does 
not follow this pattern, and it has the rms increasing with the increasing 
reflection. This appears to indicate that the high-rms of this soft state may 
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Figure 8: (a) The total rms vs. the amplification factor, showing the hard state 
(black filled squares), the intermediate state (magenta empty squares) and two 
kinds of the soft state, with low total rms (<0.22, blue filled circles), and high one 
(>0.22, red empty circles), (b) The total rms vs. reflection. The rms decreases and 
the reflection strength increases from the hard to soft state, except for the high-rms 
soft state (grey region), which has also unusally strong reflection. 
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Figure 9: The rms spectra, a(E)/F(E), for the six selected observations for which 
the energy spectra are shown in Fig. [3l 



be due to its strong reflection component. 

We then study the fractional variability as a function of the photon 
energy, rms(E) = a(E)/F(E), where u{E) is the standard deviation at the 
energy E. To calculate it, we use the method described in Gierlihski &: 
Zdziarski (2005), in which the light curves are binned at (1/256) s. Figs. 
EKa-f) present rms(E) for the six selected observations. We see rms(-E') is 
constant in the hard state, decreasing with energy in intermediate states, 
and increasing with energy in soft states. 

We can model the rms(E) by three variability patterns (Gierlihski &: 
Zdziarski 2005). The simplest one is given by a varying overall normaliza- 
tion, which can correspond to a varying accretion rate at a stable geometry 
(with an optically-thick disc truncated at a large radius), see Fig. [ToT a). 
This fits very well the hard-state rms of Fig. El^a), as shown in Fig. PTOT b). 
The second pattern (applicable to hard/intermediate states) corresponds to 
a variable irradiation by disc photons and the constant inner flow, which we 
model as varying £ s and constant £^, see Fig. HOfc), This fits well the rms in 
Fig. E]^c) (decreasing with energy), as shown in Fig. [TOTd). Then we consider 
a constant disc and a variable hot corona, which we model as varying £^ and 
constant £ s , see Fig. PTOl e). This fits very well the soft-state rms of Fig. [9^f ), 
as shown in Fig. [TUlff). 

3 Discussion and conclusions 

Considering the energy spectra of Cyg X-l, we have found good fits of the 
hybrid Comptonization (eqpair) model to all the spectra. This indicates 
possible presence of non-thermal acceleration in all of the spectral states. 
The spectral state is strongly and uniquely related to the Comptonization 
amplification factor. In the hard state, the seed disc photons are weak and 
the amplification factor is high. In the soft state, the disc photons dominate, 
the corona is relatively weak, and the amplification factor is low. We also 
see a strong reflection-index correlation, in agreement with previous work. 
The ionization of the reflecting medium increases from the hard to the soft 
state, as expected theoretically. 

Regarding the power spectra, they are well modelled by multiple 
Lorentzians. The peak frequencies of all of the Lorentzians increase from 
the hard to the soft stae. In particular, we find a very good correlation of 
£h/£ s oc fi . An additional, broad-noise (an e-folded power law), compo- 
nent is required in all states, and its importance increases with the spectral 
softness. We also find a peculiar high-rms soft state, which is characterized 
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Figure 10: (a) The theoretical spectral variability pattern due to a changing M. 
(b) The fit of this pattern to the rms spectrum of Fig. |9ja) (hard state), (c) The 
theoretical spectral variability pattern due to a changing £ s and constant £h- (d) 
The fit of this pattern to the rms spectrum of Fig. M, c ) (intermediate state), (e) 
The theoretical spectral variability pattern due to a changing £^ and constant l s . 
(f) The fit of this pattern to the rms spectrum of Fig- Eff ) (soft state). 



by unusally strong Compton reflection, and which appears not to be seen in 
other black-hole binaries. 

Regarding the rms spectra, they are difficult to interpret uniquely. Still, 
the simplest interpretation appears to be that the hard state variability is 
due to a varying accretion rate, that in the intermediate state is due to 
varying disc photon input (irradiating the hot inner flow), and that in the 
soft state, due to a variable hot corona above a stable disc. 

Similar rms variability patterns are found in transient low-mass X-ray 
binaries (Gierlihski & Zdziarski 2005). The variability patterns of Cyg X-l 
on long time scales based on the ASM/BATSE data (Zdziarski et al. 2002) 
are also fairly similar to those found here. 

Our results are in very good overall agreement with the model of a 
truncated outer disc and an overlapping hot inner flow (e.g., Done et al. 
2007). The disc moves in from the hard to the soft states, as indicated by 
the increasing spectral softness, reflection strength and the characteristic 
frequencies. 

In the hard state, the outer disc is truncated at a relatively large radius, 
and the overlap with the hot flow is weak or absent. Then, the disc moves 
in intermediate states, probably reaching the innermost stable orbit in the 
softest state. In addition, our results show how the two components of this 
accretion flow vary with the changing accretion rate. In the hard state, 
mostly the overall normalization changes, which is, presumably, due to a 
changing accretion rate. In intermediate states, the variability is driven by 
variable irradiation of the inner hot flow (with an approximately constant 
power) by soft photons from the outer disc, which inner radius varies. The 
disc moves towards the last stable orbit in the soft states and the hot flow 
becomes a hot corona above the disc. The variability is then driven by 
changing dissipation in the corona, whereas the disc remains stable. In the 
last two variability patterns, the variability is due to changing of the weaker 
spectral component, blackbody emission from the outer disc in the hard 
state, and hard X-rays from the hot plasma above the inner disc in the soft 
state. 
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